Chemical Potential Jump during Evaporation of a Quantum 

Bose Gas 



o 



- 1—1 

X 



E. A. Bedrikovalll A. V. LatyshevE 



Faculty of Physics and Mathematics, 
Moscow State Regional University, 105005, 
Moscow, Radio str., 10- A 



m 
o 

(N 

Ch ■ Abstract 



The dependence of the chemical potential jump coefficient on the evaporation 
coefficient is analyzed for the case in which the evaporating component is a 
' r~| '■ Bose gas. The concentration of the evaporating component is assumed to be 

much lower than the concentration of the carrier gas. The expression for the 
chemical potential jump is derived from the analytic solution of the problem 
for the case in which the collision frequency of molecules of the evaporating 
component is constant. 



PACS numbers: 05.20.Dd Kinetic theory; 05.30. Jp Boson systems; 71.35.Lk 
Q\ ■ Collective effects (Bose effects, phase space filling, and excitonic phase transitions) 



INTRODUCTION 

Interest in problems associated with the behavior of gases under 
conditions in which quantum properties cannot be ignored has been 



growing in recent years pQ. The behavior of mixtures of such gases is 
also of considerable interest. The "quantum" gases that are studied 
most frequently are 3 He n 4 He. It should be noted that 3 He is a Fermi 
gas, while 4 He is a Bose gas. Such a combination of different quantum- 
mechanical statistics attracts considerable attention to the mixtures of 
these gases [2]. A large number of publications are also devoted to the 
properties of solutions of these gases 0,11] • At the same time, boundary 
value problems for such mixtures have been studied insufficiently. Such 
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problems include the behavior of a mixture of quantum gases in the 
vicinity of the evaporation boundary. Let us consider the most frequently 
encountered case of a dilute mixture. Let molecular concentration n\ of 
one gas be much lower than molecular concentration n<i of the other 
gas: n\ <C ri2- We will consider the problem of evaporation of the 
first gas to the gas mixture. We will henceforth consider the problem 
in a more general formulation, taking into account the possibility of 
other applications also. Let us suppose that we have a mixture of two 
quantum gases. We consider the problem of evaporation of a Bose gas 
to this mixture, provided that its concentration is much lower than the 
concentration of the other component. 

1. FORMULATION OF THE PROBLEM AND BASIC 

EQUATIONS 

Let us consider evaporation from a plane surface to a binary gas 
mixture. We assume that concentration n\ of the evaporating component 
of the mixture is much lower than concentration n<i of the nonevaporating 
component (n\ <C Ti2] dilute mixture). It should be noted that this 
condition holds in the most important applications. 

In the semiclassical approximation, the Boltzmann equation for a 
binary gas mixture has the form [5] 

9fi dfi 

-Zr + Vi-z- = Ju + Jij, i^j,ij= 1,2. (1) 
at oyi 

Here, fa is the distribution function for the zth component of the 
mixture and Ja and are the integrals of collisions of molecules of 
the ith component with one another and with molecules of the jth 
component, respectively. 

It should be noted that J\\ ~ n\ and J12 ~ nmi (since fa ~ rij). 
Quantity e = n\jn<i is a small parameter [e <C 1), since n\ <C n<i. 
Obviously, | Ju\/\ Jn\ ~ £■ Consequently, quantity J\\ can be disregarded 
as compared to J12 in the first approximation in e. In addition, the 
action of the first component on the distribution function for the second 
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component can also be disregarded in this approximation in e. Therefore, 
under the conditions of the given problem, the distribution function 
for the second component of the gas mixture can be treated as an 
equilibrium function with mean velocity U2 = and constant tempera- 
ture T and concentration 112- 

Quantity J 12 can be approximated by a kinetic model of the Bhatna- 
gar— Gross— Krook (BGK) type PE]. Taking into account Eq. (1), we 
can write the kinetic equation for the first component in the form 

^ + vi^ = z/i(/ eg -/), (2) 



where 



eq 



f eq is the equilibrium Bose distribution function for the first component, 
k is the Boltzmann constant, V\ is the effective collision frequency for 
molecules of the first component, /i(r) is the chemical potential of the 
first component, and T is the temperature of the mixture, which is 
assumed to be constant. 

Let us choose the Cartesian system of coordinates with the center on 
the surface from which evaporation takes place. We direct the x\ axis 
along the normal to the surface. During evaporation from the surface, a 
constant gradient of the first component concentration 

dn\ 

Xl=OQ 
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exists at a large distance from the surface. 
Here 

n(xi) = J f{x h v) dil : dlt = — j^Kf — ' 

is the concentration (number density) of the first component and p is 
the momentum of molecules of the first component. 

We assume that evaporation is weak; in other words, we assume that 
the relative variation in the concentration of the first component over 
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the mean free path I of molecules is much smaller than unity: 

\G n \ <C 1, G n = — \g n \. 

n s 

Here, n s is the concentration of saturated vapor (gas) of the first 
component on the evaporation surface, which corresponds to surface 
temperature T s — T. 

Under such conditions, the problem can be linearized. Preliminarily, 
we pass to dimensionless velocity C = v/vt, dimensionless coordinate 
x = x\/It and dimensionless chemical potential a{x) = /i(x)/(kT). 

Here vt = ^ s the thermal velocity of the first component, 

(3 = m/(2kT), I = vtt is the mean free thermal path of molecules of 
first component, r = 1/Vi is the mean time between two collisions of 
molecules of the first component. 

We can now write Eq. (2) in the form 

C x ^ = f F (x,C)-f(x,C). (3) 

Here fs(x, C) is a locally equilibrium Bose distribution function, 

f B (x, C) = \ - 1 + exp(C 2 - a(x)) 

We will linearize the problem relative to the absolute Bosean 

f B = f B (C, a a ) = [ - 1 + exp(C 2 - a s )] \ 

where a s is the value of the dimensional chemical potential corresponding 
to surface temperature T = T s and the concentration of the saturated 
vapor at this temperature. 

We will linearize the dimensionless chemical potential relative to 
quantity a s ; i.e., we assume that a(x) = a s + a(x). Linearizing the 
locally equilibrium Bose distribution function relative to /Jj, we obtain 

f 8 B(x,C) = f a B (C,a 8 ) + g(C,a a )a(x), 



where 

g{C, a s ) = exp(C 2 - a 8 )(-l + exp(C 2 - a s )y 2 . 



We will seek the distribution function in the form 

C) = f F (C, a a ) + g(C, a s )h(x : C x ). (4) 

Using this expression (4), we can write Eq. (3) as: 

dh 

C x — = a(x) - h(x, C x ). (5) 

We can now find the dimensionless relative deviation a(x) = a(x) — 
a s . of the chemical potential from its value at the wall using the law of 
conservation of the number of particles: 

J f(x,c) dn = J f B (x,c)dn. 

The law of conservation of the number of particles leads to the equation 
J [h(x, C) - a(x)]g(C, a s ) d s C = 0, 

which gives 

/ h(x,C x )g(C,a s )d 3 C 

a(x) = J - j . (6) 

J g(C, a a ) d 3 C 

We can easily find that 



oo 



J g(C,a s ) d 3 C = 47r J g{C, a s )C 2 dC = 2nf Q {a s ), 



where 

oo oo 



dC 



foM = / /b(C, ol 8 ) dC = — — 

J J -1 + exp(C^ - a s ) 



We transform the numerator of expression (6) as follows: 

oo 

J h(x,C x )g(C, a s )d 3 C = 7i J f B (C x ,a s )h(x,C x ) dC x 



-oo 
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Consequently, in accordance with relation (6), the deviation of the 
dimensionless chemical potential is given by 

oo 

a ( x ) = WT7 — T / fB{n,a s )h{x,n)dn. (7) 
2/o(a a ) J 

— oo 

Thus, kinetic equation (5) for the problem of evaporation, taking into 
account relation (7), has the form 

00 

fi— + h(x, fi) = — — — r / f B (/i',a s )h(x, fi')dfi', ji = C x . (8) 
ox 2/o(a s ) J 



2. FORMULATION OF THE BOUNDARY CONDITIONS 

Let us consider the boundary condition at the evaporation surface 
for molecules of the first component taking into account the effect of the 
surface by introducing evaporation coefficient q (see [HIE]), 

f(0,C) = qf s B (C,a s ) + (l-q)f B (C,a \ C x > 0, (9) 

where 



f B (C,a )= -l + exp(C 2 -a 



2 



-1 



Quantity a® can be determined from the nonpercolation conditions 
for molecules reflected from the surface without being condensed on it 
(the probability of such a process is 1 — q), 

(1-9) / C x [fUC,a )64C x ) + f(0,C)64-C x )]dn = 0, 

where 6 + (x) ) is the Heaviside function, 6+(x) = 1 for x > and 
0+(x) = for x < 0. 

Taking into account the definition of the Heaviside function, we can 
transform the nonpercolation condition to 

J C x f B (C,a )d 3 C+ J C x f(0,C)d 3 C = 0. (10) 

C x >0 C x <0 
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The boundary condition for function h(x, C x ) at the wall can be 
derived from condition (10). Substituting function (4) into (10), taking 
into account condition (9) as well as the results of the linearization 

1 1 



mC,a ) = 



— 1 + exp(C 2 — ao) —1 + exp(C 2 — a s — (ao — a s)) 
= f s B {C, a s ) + g(C, a s )(a - a s ). 



and carrying out the substitution, we obtain the following boundary 
condition at the wall: 

h(0,ti) = (l-q)(a Q -a 8 ), fi>0. (11) 

At a large distance from the surface (outside the Knudsen layer having 
a thickness on the order of the mean free path of molecules), function 
h(x, C x ) has the form 

h(x,fj) = haa(x,fj) + 0(1), C x = fJL, X ->• +OO, (12) 

where 

h a s(%, 1^) = A a + G a (yX — fl), 

G a here, G a is the gradient of the dimensionless chemical potential 
defined far away from the wall, 

da(x) 



G = f da{x) \ 
V dx J 



x=+oo 



and A a is the chemical potential jump, viz., the unknown quantity to 
be determined from the solution of the problem. 

Substituting function h as (x,/i) into definition (6) of the dimensional 
chemical potential, we obtain the asymptotic distribution of this potential: 

a as (x) = A a + G a x, x — > +00. 

It follows, hence, that A a = a as (0) = a as (0) — a s , i.e., the chemical 
potential jump is defined as the difference between the extrapolated 
value of chemical potential at the wall and its value immediately at the 
wall. 
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Quantity ao — a s can be determined from nonpercolation condition 
(10). This condition can be written in the explicit form 



C x >0 



+ J C x 

C x <0 



f s B (C,a s ) + g{C,a s ){a - a s ) d?C+ 



f s B (C,a s )+g(C,a s )h(0,C x ) d 3 C = 0. 



It should be noted that the sum of the integrals of the first and second 
terms in each square bracket is zero; therefore, we arrive at the following 
equation: 

(a -a s ) J C x g(C,a s )d i C + J C x g(C,a s )h(0,C x ) d 3 C = 0. 



c x >o 



C x <0 



The second integral, which can be evaluated over the negative half- 
space, will be calculated using the conservation law for momentum (to 
be more precise, the x component of the momentum). For this purpose, 
we replace h(0, C x ) from the first integral on the right-hand side of the 
obvious equality 



J C x g(C,a s )h(Q,C x )d 3 C = 



c x <o 



= J C x g(C,a s )h(0,C x )d 3 C- J C x g(C,a s )h(0,C x ) d 3 C 

c x >o 

by h as (0, C x ) and h(0, C x ) from the second integral by (1 — q)(ao — a s ) 
from boundary condition (11). After evaluating the required integrals, 
we arrive at the equation, which gives 

4g 4 (a s ) G a 



ao - a s = 



(13) 



Here 



oo 



9n + 2(a s ) = J C n+2 g(C,a s )dC, 



n = 0,1,2. 
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Evaluating the first two of these integrals, we obtain 

1 3 
g3{a s ) = -- ln(l - e as ), gt{a a ) = -l{a s ). 



l(a s ) = - J \n{l-e as - c2 )dC. 



Consequently, equality (13) can be written in the explicit form 

2l(a s ) Ga 
a ° as ~ ln(l - e«.) ' 9 • 

Now, boundary condition (11) is complete 



h(0,fi) = B, p>0, r A e B = -LJL . — G(X . (14) 

q m(l — e as ) 

Thus, the boundary value problem involves the obtaining a solution 
to Eq. (8), which satisfies boundary conditions (14) and (12). Quantity 
a s will be henceforth denoted as a. 



3. DIFFUSION COEFFICIENT AND THE MASS FLOW 

RATE 

Function 

fas(x, C) = f s B (C, a) + g(C, a)h as (x, C x ), 

where h as {x, fi) = A a + G a (x — //), is known as the Chapman— Enskog 
distribution function [3,[B]. 

Using this function, we can calculate diffusion coefficient D12. Diffusion 
flux i emerges due to the presence of a density gradient Vp in the gas; 
consequently, for small values of V/?, we have i = —D^Vp [TP] . In 
accordance with the formulation of the problem, this gives 

D i_x_ _ m J f as v x rffl 

12 P'x x m J(fa S )' Xl dtt' 
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Here, we replace dimensionless coordinate x by dimensional coordinate 
x\ = x/(viyf]3) and pass to integration with respect to dimensionless 
velocity. This gives 

AkT g/t{ot) 



Dl2 = 



15 x 



3mz/i f (a) 

It should be noted that this result is transformed into the classical 
one for a — > — oo [9]. Indeed, using the asymptotic form 



TT 



a 



-oo, 



we obtain the known result 



D 



12 



kT 



It should be noted that, in the approach used here, diffusion coefficient 
D\2 is treated as an empirical quantity. For example, for the 3 He - 4 He 
mixture at a temperature of 2K, we have [11] D12 ~ 10~ 3 cm 2 /c. 

Let us find the mass flow rate of the evaporating component of the 
binary gas in the direction of the By definition, the mass flow 

rate in the x direction is 

U - = ^j fv*^ N = J f d ^ 

Substituting this expression into the definition of diffusion coefficient, 
we obtain 



f as v x dtt 



N 



N 



= Ua- 



N 



asJx 



dn 



J \fas) Xl 



dn 



It can easily be seen that 



Dn = U X -N 



-1 



passing to dimensionless mass flow rate W x = y/j3U x , we obtain 

Wt 2kT tt 2 N 

D10 = 



G a v\m (2s + l)k^g2(a) 
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where 

oo 



fo(a) = J f B (C,a)dC, k T = 



ttlvt 



h ' 



where kr is the thermal wave number. 

The number density in the linear approximation is given by 

N = j fB(C,*)dn= (2s +^ h(a), 



where 

oo oo 



/>(<*) = / f B (C,a)C*dC = J _fl%_ a = 





Thus, we obtain 

_ _ W x 2kT f 2 (a) 

u n - yr ■ i-T 

G a vim g 2 {a) 

whence the mass flow rate is 

_^rn f (a) 
Wx ~ 2kT Dl2 l{a) Ga - 



4. SOLUTION OF THE PROBLEM 

We will seek the solution to Eq. (8) in the form 

x 

h v (x,fi) = exp( — )$(7/,/i), 
r\ 

oo 

f B (fJ.,a)^{rj,fj)dfjL = 1, 



2/o(«) 

-oo 

where 77 n, is the spectral parameter or the separation parameter. 

Using these two equalities, we obtain from Eq. (8) the characteristic 
equation 

(77 - /i)$(r/,/i) = rj. 
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For rj 6 (—00,00) we obtain the eigenfunctions of the characteristic 
equation, 

i)-V> jB{r),a) 

Here, the symbol Px" 1 indicates the principal value of the integral of 
x~ , is the Dirac delta function, and A(z) is the dispersion function 
of the problem, 

00 

77 /" f B (fi,a)dfi 

2/o(a) J A 4 - ?7 

—00 




Phc. 1. Dispersion function on real axis, curve 1 correspond to Bose-gases, curve 2 
correspond to Fermi-gases from [15], a = —1. 



We will seek the solution to problem (8), (12), and (14) in the form 
of an expansion 



30 



/X 
exp( — )§(<q,n)A(rj) dq, 



(16) 
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where A{rj) is an unknown function (coefficient of the continuous spect- 
rum) and A a is an unknown constant (coefficient of the discrete spect- 
rum). It will be shown below that coefficients A a and G a are connected 
by a linear relation, 

A a = C(a,q)G a , (17) 

in which function C(a,q) is the coefficient of the chemical potential 
jump. 

Let us find the coefficients of the discrete and continuous spectra. 
The substitution of eigenfunctions into expansion (16) and of expansion 
(16) into boundary condition (14) leads to a singular integral equation 
with the Cauchy kernel: 



oo 

B = A a -G a »+ fUM^l^ + im. A^AOO = 0, „ > 0, 

o 



We introduce an auxiliary function, 



l 



M(z) = J 



r]A(r]) drj 



1} — z 
o 

with upper and lower boundary values on the positive semiaxis related 
by the Sokhotsky formulas 

M + (/i) - M~(n) = 2miiA(n), \i > 0, (18) 

l 

r]A(r]) drj 



i[M+(^) + M-(fi)} = MM, MM = J 



T]-fi 


and the last integral in which is treated as a singular integral in the 
sense of the Cauchy principal value. 

Using the upper and lower boundary values of functions X(z) and 
M(z) on the cut (0, oo), we can reduce the singular equation to the 
Riemann boundary value problem: 

X + ( f j)[M + ( f j) + A a -G af A-B] = 



14 



= A-(/z)[M-(/z) + A a - G a fi - B], m > 0. 

Let us consider the corresponding homogeneous boundary value prob- 
lem 



x-(ji) \-W 

Its solution |d2] has the form 



X(z) = - exp V(z), V(z) = - [ 

Z 7T J 



1 f C(r, a) dr 



r — z 



where 

c(T , a) . arcctg ?^m_.. 

7Trf B (r, a) 

Using the homogeneous problem, we can reduce the inhomogeneous 
problem to determining the analytic function from its zero-point jump 
at the cut: 

X+( A i)[M + ( /U ) + A a - G af i -B} = 
= X-(ji)[M-(ji) + A a - G af i -B} : fi> 0. 
The solution to this problem has the form 

M(z) = -Aa + G a z + B- 

A{z) 

Using the condition M(oo) = 0, we obtain the jump of the dimen- 
sionless chemical potential, 

oo 

A a = Vi{a)G a + B, where Vi(a) = ~ J C(r, a) dr. (19) 

o 

Coefficient A(tj) of the continuous spectrum can be determined by 
substituting the solution to the boundary value problem into equality 
(18): 



A{rj) = ° 



2irir] 



1 1 



-Gr 



lX+(rj) X-(n)\ 7rr]X(rj) a ' 
Thus, the unknown coefficients in expansion (16) are determined. 
This means that the distribution func-tion for the evaporating component 
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is constructed completely and can be written in the form 
h(x, ji) 



oo 



— = C (a, q) + x - M + - J exp { - - ) — sm 



Using the contour integration methods for boundary x = 0, we obtain 



MM) = ^ _ C0S C<» 



^1 — 6 + (fi)^j, — oo < fi < +oo. 



This expression shows that for fi > 0, the distribution function for 
molecules reflected from the wall exactly satisfies boundary condition 
(14). 



5. CHEMICAL POTENTIAL JUMP AND PROFILE 

Let us compare equalities (17) and (19). Using the second equality in 
(14), we obtain the coefficient of the dimensionless chemical potential: 

C( a , q ) = V l{a )-^J^ y (20) 

Figure 2 shows the dependence of coefficient K(a, q) on the evaporati- 
on coefficient q; curves 1 and 3 correspond to values of dimensionless che- 
mical potential a = —0.5, and —3, respectively. The curve 2 correspods 
to case of Fermi-gases at a = —0.5. 

Figure 3 shows the dependence of coefficient K(a,q) on the dimen- 
sionless chemical potential; curves 1, 2 and 3 correspond to the values 
of evaporation coefficient q = 0.5, 0.3 and 0.2, respectively. The curve 2 
correspods to case of Fermi-gases at q = 0.5. 

Equality (17) describes the jump of the dimensionless chemical poten- 
tial. Passing in this equality to dimensional quantities, we obtain the 
chemical potential jump in a quantum Bose gas: 

^(0) = C(a,?)/(^) 

V ax\ J xi=+oo 

where / = \j{y\^f]5) is the mean free path of molecules. 
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We express collision frequency v\ in accordance with relation (15). 
Then the chemical potential jump (in terms of dimensional quantities) 
can be calculated by the formula 



m n fdfi(xi] 

where K(a, q) is the coefficient of the chemical potential jump, defined 

as 



K(a,q) 



2 /o(«)r T// , l-q 21(a) 



VAa 



1(a) L v 1 q ln(l - e a ). 
Note that in the limit of high temperatures, when quantum properties 
of the gas can be ignored, result (20) is transformed into the following 
known result for the concentration jump in a classical gas [12]: 

C(-oo,q) = 1.0162+ v^F—^. 

The chemical potential distribution in half-space x ^ (known as 
the chemical potential profile) is defined by equality (7). Substituting 
the expansion of distribution (known as the chemical potential profile) 
is defined by equality (7). Substituting the expansion of distribution 
function (16) into (7), we find that the chemical potential profile can be 
constructed using the formula 



a (x 



oo 

= C(a,q) + x+^ J exp(-^- V(?7))smC(77). (21) 



G a tt J V r\ 



6. CONCENTRATION JUMP AND PROFILE 



Concentration profile N(x) of the gas in half-space x ^ is given by 
the equality 

N(x) = N s + n(x), 

where 

N s = J f F (C,a)dQ, n(x) = J h(x,C x )g(C, a) dQ. 

We can easily find that N s = N$l(a), where 

2tt(2s+ l)m 3 



27ThY(VP) 3 ' 
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and the deviation of the concentration from concentration n(x) of satu- 
rated vapor is given by 

n(x) = P N (x)G N , (22) 

where 

i- oo 

Pn(x) = 1 



C(a, q) + x H — / exp ( Vil)) sm ((l) 

7T J \ rj / 





g 2 {a) 



where Pn(x) is the coefficient of the concentration profile, Gn is the 
quantity of the concentration gradient (with respect to the dimensionless 
coordinate), 

UN ~ 



dx / x=+oo 

For the constant gas concentration, the chemical potential increases 
upon cooling [H], so that a — > +oo for T — > 0, while for T — > oo, 
a — > — oo. 

Consequently, the curves plotted for small values of a correspond to 
a lower effective gas temperature. 

It follows from Eqs. (21) and (22) that the profiles of chemical potential 
and concentration are proportional to within the corresponding gradients. 

In addition, it can easily be seen that the gradients of chemical 
potential and concentration are proportional to each other with a factor 
Nq\ Gn = NoG a . 

The latter equality indicates that the chemical potential gradient is 
the gradient of the logarithm of concentration: 

'dlniV- 



Grv — 



x dx 



x=+oo 



We can easily find that the relation between the deviations of con- 
centration and chemical potential from the corresponding values at the 
wall is 

?P- = g 2 (a)a( x ). (23) 

iV 

Using this relation, we can find the concentration jump for x = 0, 
n as (0) = N g 2 (a)a as (0) = C(a, q)g 2 (a)G N = C N (a, q)G N , (24) 
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where Cat (a, q) = C(a, q)g2{ct) is the coefficient of the concentration 
jump. 

The behavior of the concentration jump coefficient from relation (24) 
is depicted in Figs. 4 and 5. 

Figure 4 shows the dependence of the concentration jump coefficient 
of the evaporation coefficient; curves 1 and 3 correspond to dimension- 
less chemical potential a = —0.5 and a = —3, respectively. The curve 
2 corresponds to case of Fermi-gases at a = —0.5. 

Figure 5 shows the dependence of the concentration jump coefficient 
on the chemical potential a; curves 1,3, and 4 correspond to values of 
accommodation coefficients q = 0.5, 0.3, and 0.2. The curve 2 corresponds 
to case of Fermi-gases at q = 0.5. 

7. CONCLUSIONS 

The problem of evaporation of one of the binary gas components has 
been solved analytically. The evaporating component is a Bose gas. We 
studied the dependence of the chemical potential jump coefficient on the 
evaporation coefficient and on the chemical potential. On the basis of the 
analytic solution of the problem, an explicit representation is obtained 
for the distribution function, as well as the expressions for the chemical 
potential jump and the chemical potential distribution in the half-space 
(chemical potential profile). It is also shown that the concentration jump 
for the Bose gas and its distribution in the half-space are proportional 
to the chemical potential jump and distribution. On Figs. 2-5 graph 
comparison of the received results with similar results from work [TB] is 
shown . 



19 




Phc. 2. The dependence of coefficient K(a, q) on the evaporation coefficient q; curves 1 
and 3 correspond to values of dimensionless chemical potential a = —0.5, and —3, 
respectively. The curve 2 correspods to case of Fermi-gases at a = —0.5. 




Phc. 3. The dependence of coefficient K(a, q) on the dimensionless chemical potential; 
curves 1,2 and 3 correspond to the values of evaporation coefficient q = 0.5,0.3 and 0.2, 
respectively. The curve 2 correspods to case of Fermi-gases at q = 0.5. 
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Phc. 4. the dependence of the concentration jump coefficient of the evaporation 
coefficient; curves 1 and 3 correspond to dimensionless chemical potential a = —0.5 and 
a = —3, respectively. The curve 2 corresponds to case of Fermi-gases at a = —0.5. 




Phc. 5. the dependence of the concentration jump coefficient on the chemical potential a; 
curves 1,3, and 4 correspond to values of accommodation coefficients q = 0.5,0.3, and 
0.2. The curve 2 corresponds to case of Fermi-gases at q = 0.5. 
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